Exclusive Neutral Pion Electroproduction in the Deeply Virtual Regime 
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We present measurements of the ep — >■ epTr" cross section extracted at two values of four- 
momentum transfer — 1.9 GeV^ and — 2.3 GeV^ at Jefferson Lab Hall A. The kinematic 
range allows one to study the evolution of the extracted cross section as a function of and W. 
Results are confronted with Regge-inspired calculations and GPD predictions. An intepretation of 
our data within the framework of semi-inclusive deep inelastic scattering is also discussed. 

PACS numbers: 13.60.Hb, 13.60.Lc, 13.87.Fh, 14.20.Dh, 25.30.Rw 



I. INTRODUCTION 

The past decade has shown a strong evolution of 
the study of hadron structure through exclusive pro- 
cesses, allowing access to the three-dimensional structure 
of hadrons. Exclusive processes include deeply virtual 
Gompton scattering (DVGS) and deeply virtual meson 
production (DVMP). This document focuses on the lat- 
ter, and more precisely on neutral pion production. 
We present measurements of the differential cross sec- 
tion for the forward exclusive electroproduction reaction 



ep — > epTr", through virtual photoabsorption. A dia- 
gram of this process, including definitions of the kine- 
matic variables, is presented in Figure [T] Results will 
be presented for four kinematics. Two of them are de- 
fined by the same value of xbj = 0.36 and are called 
Kin2 (at = 1.9 GeV^) and Kin3 (at = 2.3 GeV^) 
The two remaining ones are defined by the same value 
of g2 = 2.1 GcV^ and are called KinX2 (at xsj = 0.40) 
and KinA3 (at xbj = 0.33). 

The behavior of the cross section will be compared to 
different models that are available to describe tt^ electro- 
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FIG. 1: Diagram of the forward tt^' electroproduction reaction 
(top), and of the dominant tt" decay mode (bottom). The 
kinematic invariants of this reaction are defined as: — 

,c.m.| \ICM\\2 



-{k - k'Y, XBi = QV{2pq), t = {q-^q' 
Q^l/xBj - 1), and t^in = ('5'-'""'' 
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^'CM\ ^j^g norms of q, q' in the pv^ final 
state center-of-mass frame. 



production, including the Rcgge model and the general- 
ized parton distribution (GPD) framework. 

Forward photo-production at asymptotically high en- 
ergies can be described by the Regge theory, which ex- 
ploits the analytic properties of the scattering amplitude 
in the limit t/s Previous analyses have applied 
Regge phenomenology to exclusive photo- and electro- 
production in the kinematic range presented here [H, Q . 
Recent computations with Regge-inspired models exist 
for our kinematics. These models include p, uj, and 
h meson exchange as well as tt^ rescattering. Among 
these, there is the t-channcl meson-exchange (TME) 
model by Laget et al. A brief description of this model 
has been given in [3], and it is described extensively in 
[^, Recent JLab Hall C experiments studying the 
dependence of charged-pion electroproduction with 
a longitudinal-transverse separation were analyzed using 
the TME formalism ^ . Another Regge-inspired compu- 
tation by Ahmad, Goldstein and Liuti 0] is available for 
our kinematics. 

In the Bjorken limit oo, and t/Q^ <^ 1 at fixed 

XBj , the scattering amplitude is dominated by the lead- 
ing order (or leading twist) amplitude of GPDs and the 
pion distribution amplitude (DA) The GPDs are 

light-cone matrix elements of non-local bilinear quark 
and gluon operators [T2m^ . unifying the elastic elec- 
troweak form factors with the forward parton distribu- 
tions of deep-inelastic lepton scattering. Cross section 
predictions within the GPD framework exist for the lon- 
gitudinal cross section (Jl [13, EH ■ With the definitions 
of the cross sections are predicted to scale as 

o'L ~ and <tt ^ . Thus at sufficiently high Q^, 
(jl will dominate over ar- Beam spin asymmetries for 
forward exclusive 7r° electroproduction have been mea- 
sured for > 1 GeV^ [l] . We performed measurements 
at two values at fixed xb in order to test these predic- 
tions of dependence. An interpretation of exclusive 
data with semi-inclusive mechanisms also exists to ex- 
plain transverse cross sections of hard exclusive charged- 
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FIG. 2: (Color online) Distribution of H{e,e'Tv'^)X events in 
the [xBi, Q'] plane, for Kin2 (xBj = 0.36, = 1.9 GeV^) 
and Kin3 (xbj = 0.36, = 2.3 GeV^). Events for KinX2 
{xBj = 0.40, = 2.1 GeV^) and KinX3 (a;Bj = 0.33, = 
2.1 GeV^) are bounded by the two horizontal lines. 



pion electroproduction [l5|. 

In the second section details of the experiment are pre- 
sented, while the third section is devoted to the calibra- 
tion of the calorimeter. The formalism of 7r° electropro- 
duction by Drechsel and Tiator [l^l is presented in the 
fourth section with a special emphasis on the expressions 
for the hadronic tensors. The fifth section is devoted 
to the extraction of the cross sections and the sixth and 
seventh sections to the radiative corrections and the eval- 
uation of the systematic errors. Finally, our results are 
presented in Sec. IVIIIl with a discussion and conclusions 
in Sees. IIXI and 1x1 respectively. 



II. EXPERIMENT 

The present data were acquired as part of Jefferson 
Lab HaU A experiment EOO-110 [l^. Additional details 
about the experimental co nfig uration, calibrations, and 
analysis can be found in jlSl. [Tgj. This paper reports 
on the analysis of the triple coincidence H{e, e'jj)X 
events. A 5.75 GeV electron beam was incident on a 
15 cm liquid hydrogen target, for a typical luminosity of 
lO'^^ cm^'^ s^^. Electrons were detected in a high resolu- 
tion spectrometer (HRS). Photons were detected in a 132 
element PbF2 calorimeter, each of the elements measur- 
ing 3 X 3cm^ X 20X0. The high resolution allows one to 
accurately define (1) the virtual photon, having the kine- 
matics centered at a fixed xbj = 0.36 and two values of 
Q2 ^ 1.9 and 2.3 GeV^, as shown in Figurc[2]and (2) the 
real photon momentum unit vector, thanks to the vertex 
resolution of the HRS, and the position resolution of the 
electromagnetic calorimeter. The validation threshold for 
the data acquisition trigger was set to about 1 GeV for 
each photon cluster. For the exclusive 7r° — >■ 77 events, 
the minimum distance between the centroids of the two 
clusters that guarantees separation is about 10 cm. This 
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is achieved by the minimal opening angle « 2m,r/£'7r and 
the distance from the center of the target to the calorime- 
ter front face L = 110 cm. The achieved coincidence re- 
solving time between the scattered electron and either 
photon cluster is 0.6 ns, rms. 

Figure [2] shows the distribution of i?(e,e'7r°)X events 
in the [xbj , Q^] plane, for missing mass squared = 
{q -\- p — q'Y < 1-15 GeV^. The analysis relies only on 
two specific qualities of the experiment: 

i Thanks to the resolution of the spectrometer and 
the calorimeter, one can use the missing-mass 
squared to ensure exclusivity. The exclusive sam- 
ple is selected by putting a cut on the missing-mass 
squared at the proton plus the pion mass squared. 

ii For exclusive events, the reconstruction of the in- 
variant momentum transfer t and tmin relies on the 
positions of the reconstructed photons, of which 
the resolution is better than that of the energy. 
From this, a resolution in t better than that in 
the energy is obtained. All data are presented as 
a function of imin — t, which is directly linked to 
the angle of the pion production relative to the vir- 
tual photon direction in the center of mass 9'^™ -. 

- t = 2g<= '°-g'^^=f(l - cos6'^"^ ). 

In the ep — > e'7172^ reaction, there are six four- 
vectors, equivalent to 24 independent kinematic vari- 
ables. The measured four-vectors k, p, and k' , and four- 
momentum conservation, reduce the number of indepen- 
dent variables to eight. The measurement of the two 
directional vectors ^(71) = gl/gi and ^(72) = 92/(72 
from the target vertex (reconstructed by the HRS) to 
the two cluster positions in the calorimeter provides four 
more kinematic constraints. Finally, the hypothesis that 
the observed calorimeter showers are due to photons 
(rrigj = = 0) provides two more kinematic con- 

straints. The remaining two unknowns, which we express 
as m^^ = -I- 92)^ and M|., are determined by the pre- 
vious constraints plus the energy of the two photons. Fig- 
ure [3] displays the distribution of the H{e, e'jj)X events 
in the [M^,, m^^] plane, for KinS. The upper left panel 
of this figure shows a clear correlation between the two 
variables in the exclusive region (M|- ~ Mp). This is 
a consequence of resolution fluctuations in the energies 
El and E2 of the two photons issued from a tt", which 
correlate fluctuations in M"^ and m-y^ . The missing mass 
in the right-hand panels is obtained by an empirical ad- 
justment: 

Ml] = Mil + Cim^y - m,r), (1) 

I coiT ^ I raw \ I ! " / ' \ / 

with C 13 GeV. 

This transformation produces a noticeable improve- 
ment in the distribution (lower right panel of Figure 
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FIG. 3: [(a),(b)] Distributions of H{e,e'^^)X events within 
cuts in tlie [Mx, m^-y] plane for KinS. (a) Raw distribution 
showing a clear correlation between these two variables, (b) 
The same distribution after a rotation around (M^, m^o) to 
improve the Mx resolution, [(c), (d)] Projections on the Mx 
axis of the [M%, ] distributions shown, respectively, in 
(a) and (b). The lower right panel shows that the resolution 
is indeed improved by the rotation. 



III. CALIBRATION 

Wc performed clastic H{e,e'^^i^ Phrs) calibrations at 
the beginning, middle, and end of the experiment [20j . 
The calorimeter was retracted to a position at 5.5 m 
from the target, in order to optimize the electron cov- 
erage in the calorimeter with the proton acceptance of 
the HRS. These data were used for the block calibration. 
After calibration the calorimeter energy resolution was 
observed to be 2.4% at 4.2 GeV with a position resolu- 
tion of 2 mm at 110 cm from the target. The elastic 
data also provided a consistency check on the efficiency 
of the detectors and all associated electronics from the 
observation that the elastic cross section agreed with the 
Kelly form- factor parametrization [2l| at the 1.1% level. 
During the experiment, the light output from the PbF2 
blocks decreased by up to 20%, strongly correlated with 
the distance of the blocks from the beam line. We at- 
tribute this to radiation damage of the blocks. In ad- 
dition, seven blocks, at random positions, showed much 
higher radiation damage. One explanation could be a 
poorer crystal quality of those crystals. We adjusted the 
calibration of each block, assuming an independent linear 
dose versus attenuation curve. In addition to radiation 
damage, each crystal received a pileup of low-energy pho- 
tons in random coincidence, resulting in a degradation of 
the energy resolution, and in a shift in the calibration as 
a function of its distance to the beam line. This effect 
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FIG. 4: Projection on the calorimeter of the virtual photons 
7* within cuts for Kin3. Also shown is the block relabel- 
ing used for the calorimeter calibration described in the text. 
The calorimeter is viewed from the rear, with the downstream 
beam passing to the right. 



We recall here that each event involves two blocks. The 
reconstructed missing mass of one block is then influ- 
enced by contributions from all other blocks. Because of 
this, several iterations are necessary. Then, the missing- 
mass distribution of each block for simulated events is 
adjusted to get the same missing-mass position and res- 
olution as the experimental missing-mass distribution. 
The missing-mass cut applied to ensure exclusivity is the 
same for simulation and data, and if the resolution is 
better for simulation, applying such a cut will remove 
more experimental events than simulation particularly 
near the beam where the noise degrades the experimen- 
tal resolution. This gives a spurious contribution to the 
cos (pTT term which has to be removed by smearing the 
simulation resolution. To this purpose, the momentum 
of each event i at the nth iteration contributing to the 
distribution of the block /i is changed from (<z7i)n-i 
to {qji)l, with a sampling from a Gaussian distribution: 



was taken into account through successive steps: 

i For each block the position of the reconstructed 
missing-mass squared peak was centered at 
through an energy calibration of the experimental 
data. 

ii A GEANT simulation generated a sharper resolu- 
tion in missing mass than the experimental data for 
each calorimeter block. For each block the energy 
of the simulation was calibrated together with a si- 
multaneous energy smearing, in order to center the 
reconstructed missing- mass peak position at M^, 
and to equate the resolution of the simulation to 
that of the experimental data. 

These calibrations are explained in the following para- 
graphs. 

We consider only the 90 blocks of the inner calorime- 
ter (see Figure m for the labeling), indexed by ^. We 
will assume that the energy of the photon is driven by 
the block where the shower makes the largest energy 
deposit. The 90 distributions of missing- mass squared 

(M2 ); = {k+P-k'-q^-q,)l = [Ely-{Px'y arebuih 
with all events i involving block fi. Note that for each 
event i, the reconstructed missing-mass squared appears 
in two distributions. To compare these distributions, two 
estimators are constructed: the mean {M\)^ and the 
sigma (7^1 of a Gaussian fitted to these distributions, over 
a limited range (0.62 GeV^ < (M|.)^ < 1.09 GeV^). The 
calorimeter is calibrated using 



AM|. 



-2A<z„ Ex 



Px-qj, 
\<i,\ 



(2) 



with AMI- 



((M|.)^) - M^. Neglecting the Px term 
compared to Ex between the parentheses, we obtain an 
energy correction: 



2{Exy 



(3) 



Gauss (g^) 



V2 



(4) 



where 



Aa^ = 0, 



1 2. 

simu ' 



(5) 

simu 

(6) 

and {q^)n is given by equation ([3]), except we put AM^ = 
{{{M'^)f})simu - {{M^)f,)data)/2 in this case. The factor 
2 in the denominator of AM^ is used to ensure a smooth 
convergence. 



TABLE I: Mean deviation and resolution width of the tt — > 
77 reconstruction of the data and simulation. Events are 
selected by Mx < 1-15 GeV'^ and calorimeter threshold 
-Ethr = 1.0 GeV. 





{m — m^o) (GeV) 


V((m-m,o)2) (GeV) 


Kin3 


Data 


-0.00081 


0.0088 


Simulation 


+0.00072 


0.0089 


Kin2 


Data 


-0.00017 


0.0079 


Simulation 


+0.00191 


0.0085 



The results of these iterations are shown in Figure \5\ 
Figure inland Table |T] illustrate the quality of the final cali- 
bration adjustments. The calibration of the missing-mass 
squared was cross-checked by comparing the invariant- 
mass distribution of both photons in each event. Table |T] 
lists the mean values of these distributions with respect 
to the pion mass, and their resolution. The agreement 
of the calibration with the data is at the 1.9 MeV level, 
while the widths of these distributions agree to better 
than 1 MeV. 
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the norms of (f, q' in the center-of-mass frame. 
All these quantities are defined using the convention of 
Drechsel and Tiator [l^ : z axis along the virtual photon, 
y = (kiAkf)/ sin Be orthogonal to the leptonic plane, and 
X = y A z. 

To lowest order in the fine-structure constant a, the 
differential cross section for an electron of helicity h is 



d^^ ~ dtdfjiT, ■ 



(7) 



FIG. 5: (Color online) Different iterations of the calibration 
for Kin2. The differences between simulation and data of the 
missing-mass peak position (a) and resolution (b) are shown 
before calibration (crosses), after calibration (open circles), 
and at a random iteration during calibration (asterisks). 




FIG. 6: (Color online) (a) Raw H{e, e''K^)X missing-mass dis- 
tribution for Kin3 (solid histogram) compared to the simula- 
tion (dashed histogram), and the difference between the two 
(dotted histogram), (b) H{e,e ■n^)X missing-mass distribu- 
tion at different values for the calorimeter threshold, corrected 
with a factor 1/(1 — 2(i5tiir/lj'7r|))- This correction adds to 
the distribution all tt'^ events missed because of the threshold 
value. 



IV. CROSS-SECTION ANALYSIS 

In order to extract the differential cross section, it is 
advantageous to incorporate all model-independent kine- 
matic dependences of the differential cross section into 
the experimental simulation. To this end, we express the 
differential cross section in terms of structure functions as 
described in the paper of Drechsel and Tiator [l^ directly 
related to bilinear combinations of the Chew-Goldberger- 
Low-Nambu (CGLN) helicity amplitudes [l^l- We define 
the differential phase-space elements d^^^ = dQ^dxBjd(t)e 
and = d^^ed[tmin — t]d(t>TT and the equivalent real 



(8) 



with k^ = (W^ - M^)/2Mp and k and k' the energies 
of the incident and scattered electron, respectively. The 
virtual photo- absorption cross section is expanded as 



1 



d'^av{h) _ 



{Rt + £lRl + eRrr cos2(^t 



+ \/2e^JT+l)RTL cos 0^ 

+ h^2eL{l- e)RTL' sin (j>^} 



(9) 



g" "'- = Iffl X AIp/W is the cm. virtual photon 
three-momentum, e = + 2{q^ / Q'^) ia.T? 9e/2] is the 
degree of linear polarization of the virtual photons, and 
e^/e = 4Mpx|j/(5^- The response functions are defined 
as functions of the usual hadronic tensor W^'^: 



Rt 



cos (P^Rtl = -R-eW:rz, 



sin <t)-^RTL' 

C0s2(f>T^RTT 



-ImW, 



yy 



(10) 

(11) 

(12) 
(13) 

(14) 



The interference terms Rtl and Rtl' have a leading 
sin SJr'™' dependence, and the linear polarization inter- 
ference term Rtt has a leading sin^ 6*5^™' dependence. 
For this reason, we define reduced structure functions r\, 
which remove this phase-space dependence, which are di- 
rectly related to bilinear combinations of the CGLN he- 
licity amplitudes Fi p^ : 




1 [ Rtl 
sine--- I Rti^, 



frprp 



Rtt 



sm 



rL = Rl, 



(15) 

(16) 
(17) 
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(18) 



Since our kinematics cover a wide range in xbj as well 
as in Q^, we also have to include the and the W depen- 
dence of the hadronic tensor [Wxx + Wyy)/2 + €lWzz = 
fT + ^LfL- We perform a preliminary extraction of the 
cross section on the kinematic points Kin2 and Kin3 (re- 
spectively KinX2 and KinX3) to get an estimate of the 
(respectively, W) dependence of the hadronic tensor. 
The extracted and W dependences are then intro- 
duced explicitly in the formalism to perform a second 
"definitive" extraction. The dependence is modeled in 
the form (Q^)" and . With the first iteration, the 
cross sections changed by 3%, but with a second itera- 
tion the cross sections changed by only 0.3%. 

The results will be presented as four separated cross 
sections following the usual decomposition found in the 
literature: 



dtd(j>Tr 



1 

2^ 



dax 

~dr 



d(7L 
'dt 



+ V2eL(l + e)^cos( 



dt 



e — - — cos ZcpTr 



dt 



+ hy/2eLil 



— ddTL' . 

e ) — ; — sm ( 
' dt 



(19) 



V. EXTRACTION 

We define a compact notation that summarizes Eq. 
in the form 



d^a 
d5$ 



A 



(20) 



with J^\{xy) containing all the kinematic dependence, 
A G {T+slL, TL, tt, TL'} and Xy summarizing all vari- 
ables k, Q^, XBj, W, t, considered at the vertex. T + e^L 
refiects the fact that we used only one incident energy 
and consequently, we were not able to disentangle dax 
and daL- This notation will be convenient to use for the 
presentation of the extraction process. 

The experimental data used for the analysis have the 
kinematical coverage shown in Figure [5] The analysis 
includes a complete simulation of the resolution and ac- 
ceptance of the HRS, the external and internal radia- 
tive effects on the incident and scattered electron, and a 
GEANT based simulation of the acceptance and response 
of the PbF2 array. Simulation events are generated uni- 
formly in the target vertex v along the beam line, and 
uniformly in a phase space A^$. This results in well de- 
fined values of ^^i-™' in each bin. The At bins are the 
same in the generation and experimental phase spaces, 
but resolution and radiative effects can cause the migra- 
tion of events from one bin to one of its neighbors (Figure 
[7]). Rather than extracting average cross sections in the 




tmin-t(GeVO 



7: 



(Color online) Raw H{e,e'^^)X distribution in the 
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FIG 

[tmin — t, Q''] plane with cuts for Kin3. The vertical lines 
delimit the bins we chose in fmin — t for our analysis. Su- 
perimposed is the (imin — t) resolution for each alternate bin 
showing that each bin is larger than the resolution. 



experimental bins, we use the simulation and the theo- 
retical form of Eq. ((20|) to directly extract differential 
cross sections from the experimental yields. 

We divide the acceptance into 24 equal bins in € 
[0, 2tt] and 8 bins in tnun -te[0, 0.3] GeV^ for both the 
helicity dependent and independent parts of the cross 
section. A bin jd in the kinematic variables recon- 
structed by the detector is defined by the limits (p-j^r G 

[</'(jd),</>(jd) + A(^(jd)], (tmin-0 e [{tmin-t){jd),{tniin- 

t){jd) + A(t,nin - t){jd)], etc. The statistics AN{jd) in 
a bin jd are determined by the physical cross section at 
the vertex convoluted with the detector response: 



ANUd) 



Xy)rA 



Cu / dxd I dxyTZ{xd, Xy) J"a( 



where Xy summarizes the reaction vertex variables, 
Xd summarizes the reaction vertex variables as recon- 
structed in the detector, Axd summarizes the range of 
integration for bin jd, Axy summarizes the range of in- 
tegration for all bins jy, Cu is the integrated luminosity, 
and TZ{xd, Xy) is the probability distribution for an event 
originating at the vertex with kinematics Xy to be recon- 
structed by the detector with vertex kinematics Xd- This 
expresses the effects of detector resolution, internal and 
external radiation, detector efficiency, and anything else 
that could migrate events from vertex kinematics Xy to 
the detector kinematics Xd- For the analysis and simula- 
tion, the integral is split into a sum over the bins Axy in 
the kinematic variables at the reaction vertex: 



AN{jd) 



dxd^ 



Ax 



dxyTZ{xd, Xy) y TA{xy)rA 



A 



Because the functions J-\(xy) contain the main part of 
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TABLE II: Cuts applied in the primary extraction, r is the 
value of the so-called r function. The r function defines the 
distance of the particle from the acceptance bound, and is 
jositive (negative) if the particle is in (out of) the acceptance 
23|. The M\ and 5thr optimizations are presented in Table 



Spectrometer cuts 
—6.0 cm < V < +7.5 cm 

la^HRS planci < 3.5 Cm 

(Horizontal collimator) 

12/HRS plane I < 7.0 Cm 

(Vertical collimator) 

Ifc' -PhrsI/Phrs < 4.5% 
r > +0.005 m 



Calorimeter Cuts 
— 15.0 cm < Xcaio < +12.0 cm 
Ij/caiol < 15.0 cm 



Physics Cuts 
105 MeV < 



< 165 MeV 



the dependence on the variables at the vertex, the quan- 
tity ta in a bin Axi, will be assimilated to its average 
kft^x^ = ?'t),A in this bin. Then, the last equation can be 
summarized in a vector notation: 



(21) 



with 



J Axd J Aa:„i=bin 3„ 

(22) 

We then replace the integration by a summation over the 
simulated events i: 



«e{j„,jd} 



(23) 



where the sum is over events originating in vertex bin 
and reconstructed in bin jd ■ Ngcn is the number of events 
generated in the simulation and is the total phase- 
space factor. The matrices Kj-^ are constructed from 
simulation events, summed over all events within cuts. 
We define Nd ^ N+ + TV" with N+ (iV") the number 
of counts within cuts with positive (negative) electron 
helicity. The cuts are the same for simulation and data 
(Table HI)) . The cuts and the corrections are summarized 
in Tables |TT] and IIIIl respectively. 

A is built, assuming that the statistical error on 
the simulation is much smaller than the statistical error 
of the data: 



E 



Nd 



(24) 



The minimization of with respect to the unknown 
quantities rj^ _a results in a linear system from which the 
rj^^A are extracted. To be fully consistent, one of the two 
quantities in the numerator has to be corrected for some 
instrumental systematic effects (Table HlI)) . Note that all 
vertex bins populate experimental bins, but the detector 
bin at the largest experimental bin in (tmin — t) can re- 
ceive contributions from larger values of (tmin — t), not 
generated in the simulation. Hence, although we extract 
an Tj^ _A value for the last bin, we do not include it in our 
results, its role is only to populate the lower (tmin — t) 
bins. 



TABLE III: Correction factors applied in the data analysis. 
The radiative correction factor is the combination of the vir- 
tual radiative correction factors (vertex renormalization and 
vacuum polarization) and the cut-off independent real radia- 
tion effects (Sec. IVII) . 



Correction 


Kin3 


Kin2 


Multitracks in HRS 


1.079 


1.099 


Triple cluster in calorimeter 


1.035 


1.020 


Radiative correction 


0.91 ± 0.02 


0.91 ± 0.02 



The average values of the kinematic variables Q^, e, 
xbj, W, t, ^inini e^c., iu a bin at the vertex are 



(25) 



Because the t^^ a are by construction constant over the 
bin Axt, and the integrals of J^tl, -^tt, and Ttl' cancel 
when integrating over ^^r, we can write 



GAa:, 



XyK, 



E 



K 



(26) 



These values are summarized in Table IVlIl for quantities 
independent of the (tmin — t) bin and in Table IVIIII for 
quantities depending on the (tmin — t) bin. 

Finally, the cross sections at the point Xj^, in a bin j„ 
are obtained by 



d(TA 
dt 



Ta{xj 



(27) 



The results are displayed in Tables HXl and 1x1 The first 
table shows the results for the two kinematics Kin2 and 
Kin3, which cover the full kinematic range of the ex- 
periment, resulting in two domains of different Q^, at 
constant xbj ■ The second table shows the results for the 
two kinematics KinX2 and KinX3, which only cover the 
domain between the two horizontal lines in Figure [2J in 

order to have two domains of different xbj at constant 
r)2 
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VI. RADIATIVE CORRECTIONS 



VII. SYSTEMATIC ERRORS 



The external radiative effects on tfie incident electron, 
and internal real radiative effects at the vertex are treated 
in the equivalent radiator approximation (23 . [25l | . Prera- 
diation is modeled by generating an event-by-event en- 
ergy loss AE'in of the incident electron (Eq) following a 
distribution (fe ~ 4/3): 



liniEo, ^Ein, tin) — 



btin + Ss/2 



AEu- 



with 



2a 

ds = — 
n 



En 



1 



bti„+5s/2 



(28) 



(29) 



where tin is the event-by-event target thickness (in radi- 
ation lengths) traversed by the electron before the scat- 
tering vertex. The Schwinger term Ss models the inter- 
nal pre-radiation. The scattered energy at the vertex is 
E[, = Eq — AEin — Q'^/ {2MpXBi)- Internal post-radiation 
is modeled by a similar distribution in the post-radiated 
energy AiJout: 



Ss/2 

Ai?niit, 



AEq 
EL 



Ss/2 



(30) 



These radiative effects are treated within the peaking ap- 
proximation. External post-radiation by the scattering 
electron is modeled with the GEANt3 simulation. Kine- 
matic shifts {e.g., in either the norm and direction of q) 
from external and internal radiations are fully included in 
the simulation and thereby unfolded from the extracted 
cross sections. 

In addition to these radiative effects incorporated into 
our Monte Carlo, we correct the data for internal virtual 
radiation (vacuum polarization and vertex renormaliza- 
tion effects) as well as the cut-off independent effect of 
unresolvable soft real radiation. These contributions are 
calculated by the following terms, respectively f26|: 



-'vacuum 
^vertex 



2a 
3^ 

a 

TT 



Oroal.O — — 
TT 



2 \E' 



Sp cos^ ^ 



2 

ll 2 















ml) 


3 



(31) 



where Sp(cos2 6*0/2) is the Spence function. After an 
approximate resummation, the correction we apply to the 
raw counts (to obtain the equivalent Born approximation 
cross section) is 



radcorr 



(1 - 5v 



.I2f 



(32) 



The numerical values for our kinematics arc tabulated in 

Table nn 



Two classes of inclusive hadronic electroproduction 
channels compete with the exclusive H{e, e'Tr^)p reaction: 
the H{e, e'TT'^)NTT, Ntttt, ... channels, with a threshold at 
= {Mp + m^f = 1.15 GeV^ and the iJ(e, e'7r°)7p 
channel. The first class includes N* and non-resonant 
Nix production in the final state, and diffractive p"*" — > 
TT+TT*^ production via the — > ep^n reaction. All these 
channels can be observed in a missing-mass squared dis- 
tribution (Figurc[B]). The -ff (e, e'7r°)7p channel originates 
from the diffractive ep — >■ epw reaction, with a 8.5% 
branching-ratio decay channel (27j . In our acceptance, 
the (e,e'7rO) missing-mass squared threshold for exclu- 
sive ijj electroproduction is 1.0 GeV^, thus slightly lower 
than the N-k threshold of 1.15 GeV^. However, based on 
e-p — > epo; measurements performed by [2^ , the expected 
background of W7r°7 events for M\ < 1.15 GeV^ is less 
than 1% of the exclusive i?(e, e'7r°)p yield in all tnUn — t 
bins. 

The systematic errors in the extraction method arc due 
to the cut on the missing-mass squared M\ and on the 
calorimeter threshold i?thr. The stability of the results 
is checked by varying each cut in turn. The variation in 
the estimator 



(33) 



bin=0 



is used to quantify the systematic errors. 



i For the exclusivity (M^) cut, we consider the sta- 
bility interval from 0.9 to 1.10 GeV^ in the M\ 
cut. At the high end we expect the cross section to 
have contributions from inelastic final states (Fig- 
ure (HI). At the low end, we arc removing roughly 
half of the statistics, and wc become progressively 
more sensitive to the experimental line shape. The 
stability of the exclusivity cut {e.g. for Kin3) is 
plotted in Figures [H] The cuts and variation are 
listed in Tables [TVl and fVl In each case, this study 
is performed with E'thr fixed at 1.0 GeV. 

ii For the calorimeter threshold -Bthr, the stability of 
R is expected when the software threshold is fixed 
above the hardware threshold. Above the hard- 
ware threshold, the cut is directly correlated with 
the 7r° — > 77 decay phase space, and the num- 
ber of events decreases linearly with i?thr. This 
comes from the isotropic decay of the pion, lead- 
ing to a flat energy distribution of each decay pho- 
ton. Figure [9] shows for Kin2, the quantity R along 
with the raw number of counts. The stability is 
indeed no longer observed when the statistics are 
not linear with the threshold, meaning the hard- 
ware threshold competes with the analysis thresh- 
old. The same behavior is shown for Kin3. For 
both kinematics, the systematic error coming from 
the calorimeter threshold is evaluated as ±1%. 
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c 1.45- 



+ 

1.4- 



1.35 — 



^ , I I I , I , I , I , I L 

0.9 0.95 1 1.05 1.1 1.15 1.2 

M| cut (GeV^) 



TABLE IV: Values of the cut and Bthr defining the global 
cross-section stability domain. Minimum and maximum are 
the bounds of this domain, and optimum is the cut value set 
in the middle of the stability interval. 



Variable 


Minimum 


Optimum 


Maximum 


KinS/KinXS 


Ml- cut (GeV^) 


0.90 


1.00 


1.10 


Etbr (GeV) 


1.20 


1.275 


1.35 


Kin2/KinJf2 


Mi- cut (GeV^) 


0.90 


1.00 


1.10 


Etbr (GeV) 


1.00 


1.075 


1.15 



FIG. 8: (Color online) Total cross section integrated over 
imin — t and (t>n, for Kin3, as a function of the M|- cut. The 
vertical lines indicate, from left to right, the minimal, optimal, 
and maximal Mx cut values of the stability domain. 



c 40000 

> 

a 

o 

u 20000 

.o 

E 











■ 

: 1 
- 1 
: i 






(a) 
















(b) 



0.9 0.95 1 1.05 1.1 1.15 
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FIG. 9: (Color online) (a) Total cross section integrated over 
imin — t and 4)Tr, for Kin2, as a function of Ethv The vertical 
lines indicate, from left to right, the minimal, optimal, and 
maximal Ethi values of the stability domain (see Table IIVI 
for KinS values), (b) The number of events as a function 
of EthT- The stability domain for Ethi shows the statistics 
linearly decreasing with Ethi- 



TABLE V: Experimental systematic errors. The first "To- 
tal quadratic" row shows the quadratic sum of all experimen- 
tal helicity-independent systematic errors. The second "Total 
quadratic" row shows the quadratic sum of all experimental 
systematic errors including helicity-dependent effects. 





Kin3 


Kin2 




KinX3 (%) KinX2 (%) 


Exclusivity cut 


1.5 


3.0 


Calorimeter threshold 


1.0 




HRS acceptance 


2.2 




Radiative corrections 


1.5 




Target length 


0.5 




Hadronic tensor integration 


0.3 




Multitracks corrections 


0.1 




3 clusters corrections 


0.1 




Luminosity 


0.1 




Dead time 


0.1 




Particle identification 


0.1 




Total quadratic 


3.3 


4.2 


Beam polarization 


2.0 


Total quadratic 


3.9 


4.6 



The optimal cut is set in the middle of the stability in- 
terval (see Figures [5] and [S]). The stability interval bounds 
and the optimal values for the cut and E'thr are listed 
in Table ITVl for both kinematics. 

The reduced structure functions ta are extracted at 
the optimal value of the cuts. For the structure functions 
implied in dependences, systematic errors are taken 
as the rms difference between the r\ computed at the 
optimum cuts and the computed at each of the four 
extremities of the stability domain. 

All instrumental sources of systematic errors are shown 
along with the analysis systematic errors in Table |Vl 
Since all sources of systematic errors are independent, we 
added them quadratically. This total systematic error is 
included in Tables |IX] and [Xj 



VIII. RESULTS 

The exclusive tt^ electroproduction cross section and, 
in particular, the t/)^ dependences of its separated com- 
ponents were extracted for Kin2, KinS, KmX2 and 
KmX3. Our statistics allowed us to achieve, for the (^tt- 
indepcndent cross section, a statistical precision of 3% 
for Kin2 and KinS, and of 5% for KinX2 and KinXS. 
This difference is due to the fact that we could use the 
full statistics for Kin2 and KinS, whereas less than half 
of the statistics were available for KmX2 and KmX3. 

Figure [TUl shows (Jt + ^lO'l and Figure [TT] shows cttl, 
cttt, and (Jtl' plotted as a function of tmin — t, both for 
Kin2 and KinS. Figure [T^ shows ut + ^lO'l and Figurc[T51 
shows fiTL, fTT, and utl'^ both for KinX2 and KinXS. 
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FIG. 10: (a) Separated H{e,e'-K°)p cross section ctt + ^lcl 
as a function of tmin — t for a^Bj = 0.36. Error bars represent 
statistical errors only, (b) Ratio of ax + ^lcfl for the two 
kinematics as a function of tmin — t. The fit of this ratio 
(dashed line) indicates the dependence of the cross section. 
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FIG. 11: axL (a), cttt (b), and cttl' (c) H{e,e'-jT'^)p cross- 
section components as a function of tmin — t for the two Q^- 
values. Kin2 is represented by the open circles and Kin3 by 
the solid circles. Error bars represent statistical errors only. 
The bands (light for Kin2 and dark for Kin3) show fits pro- 
portional to sin sin^ ^5r ™', and sin ^Ji- ™', respectively. 
Refer to Table HXl for more detailed cross-section values, with 
statistical and systematic errors. 



FIG. 12: (a) Separated H{e, e'-K°)p cross section ctt + ^lctl 
as a function of tmin — t for Q2 = 2.1 GeV^ Error bars 
represent statistical errors only, (b) Ratio of ar + tLi^L for 
the two kinematics as a function of t-mi-n — t. The fit of this 
ratio (dashed line) indicates the W dependence of the cross 
section. 





-0.2 

-0.4 


-0.2 
-0.4 — 
-0.6^ 
0.4 

0.2 



o W = 2.0 GeV 
• W = 2.3 GeV 



(b) 





0.05 



0.1 



0.15 



0.2 

t„i„-t (GeV^; 



FIG. 13: aTL (a), (jtt (b), and cttl' (c) H{e,e'-K^)p cross- 
section components as a function of tmin — t for the two xbj- 
values. KinX2 is represented by the open circles and KinX3 
by the solid circles. Error bars represent statistical errors 
only. The bands (light for KinX2 and dark for KinX3) show 
fits proportional to sin SJi-™', sin^ and sin ^Ji-™', respec- 

tively. Refer to Table |X] for more detailed cross-section values, 
with statistical and systematic errors. 



Wc performed fits proportional to sin 6'^™- for otl and 
cttl', and proportional to sin'^ 6'^ '"' for (Ttt- These fits, 
including statistical and systematic errors, are shown as 
bands in Figures [TT] and [T31 and in Tables IXII and IXIII 
Their reduced are below 1.05 for the Q^-dependent 
data, and below 0.75 for the xbj -dependent data. This 
confirms that the main t dependence of (Ttl,tl', and cttt 
is given by sin 9^™- and sin^ ^Ji-'™', respectively. 

The lower panel of Figure [10] (respectively. Figure fT2)) 
also shows the dependence (respectively, ccbj depen- 
dence) for the total cross section ctt + ^lCl- To investi- 
gate a or a xbj dependence, the ratio of ctt + elcl 



for the two kinematics is plotted as a function of tmin — i- 
This ratio is found to be independent of t, thus the value 
of this ratio is fitted by a constant at the zbj- and Q^- 
values for the two kinematics. 

The dependence of (Tt + clo-l in Figures [TU] and [TH 
yields the following conclusions: 

i The ratio [ctt + eLO-LlKins/lo-T + eL0-L]Kin2 is flat in 
tmin— t with a reduced of 0.33. The ratio is found 
to be 0.633 ±0.009, indicating a dependence of 
the total cross section of about 

ii The ratio [(7T + eL<yL\KinXz/WT+f-L<yL\KinX2 is also 
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flat in tmin — i with a reduced of 0.56. This ratio 
is found to be 0.660 ±0.015, indicating a W depen- 
dence of the total cross section of about l/W^^^. 

The and W dependences of the relevant quanti- 
ties [ar + clctl, ctt, and ctl, with our conventions (i.e 
Drechsel-Tiator) and VGG conventions] have been sum- 
marized in Table IVII 



Quantity 



Q dependence W dependence 



^ , , „ ^^2\-2.39±0.0S 
CTt -I- eL<7L (Q ) 

GL (Drechsel-Tiator) (q2)-0-50±o.i3 

aL (VGG) (q2)-1.50±0.08 



(T4/)-3-48±0.11 
-0.46±0.57 



{wy 

{W) 



1.28±2.52 



TABLE VI: and W dependences for the total cross sec- 
tion and the longitudional cross section with Drechsel-Tiator 
conventions and with VGG conventions. For gl, the depen- 
dences have been evaluated neglecting ctt. The and W 
dependences of ot alone {i.e. assuming gl = 0) are the same 
as the and W dependences of ctt + elctl. 
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IX. DISCUSSION 

In the domain in <min — t where we extracted cross sec- 
tions, the r\ values from Eqs. ([T5|) and are constant 
within statistics, as evidenced by the fits in Figures [TT] 
andm 

The data we extracted (see the previous section) yield 
two conclusions with regard to the available models: 

i The i-channel meson-exchange model of Laget 
(Figure [T3]) is able to describe ar + ^Lf^L and cttl', 
but neither cttl nor cttt 0- 

ii the dependence of the cross section (Figure [TOl 
and Table IVip demonstrates that we are far from 
the QCD leading twist prediction of daL/dt, which 
behaves as On the other hand, it is similar 
to the dependence of the transverse cross sec- 
tion for charged pion electroproduction published 
by Hall C 0. 

Moreover, the tt'^ has no charge and no spin, so a 
direct coupling with a virtual photon is suppressed, 
which removes the pion-pole contribution to the longi- 
tudinal cross section. This suggests that the transverse 
ep — >■ ep7r° cross section is likely to dominate, and trans- 
verse ep — >■ enTT"^ cross sections have already been de- 
scribed by quark fragmentation mechanisms usually used 
to describe semi-inclusive processes. 

T. Horn et al. measured the exclusive tt^ electropro- 
duction cross section at = 1.60 and 2.45 GeV^, with 
ctt and separation Q . The t-channel meson-exchange 
model by Laget reproduces the ctl component. How- 
ever, the Ctt component does not follow the TME model 
prediction. Kaskulov et al. performed PYTHIA-JETSET 



FIG. 14: (Color online) New calculations at Kin2 (left pan- 
els) and Kin3 (right panels) of the t-channel meson- exchange 
model, including charge pion rescattering with ttN and ttA 
intermediate states [^. Dashed lines: pole contributions and 
Pomeron cut alone. Dash-dotted lines: without pA cuts. Full 
lines: pA cuts included. 



calculations using the Lund model applied to 7r+ trans- 
verse cross sections at Hall C kinematics [l^. In this 
model, the virtual photon strikes a quark, with a proba- 
bility given by the structure functions. Due to this, the 
hadronic system fragments into two jets. The jet engen- 
dered by the single quark gives a pion, and the one en- 
gendered by the remainder of the nucleon gives the final 
neutron. These calculations applied to Hall C 7r+ trans- 
verse cross sections are in excellent agreement with the 
data. This gives evidence that the tt"*" transverse cross 
section at > 1 GeV^ above the resonance region is 
described by a partonic process. This suggests that the 
present 7r° data could similarly be described by incoher- 
ent scattering on the partonic structure of the nucleon 
target. 



For these reasons, we consider our data within the con- 
text of semi-inclusive deep inelastic scattering (SIDIS). 
We can try to fit our data with a SIDIS formalism writ- 
ten by Anselmino et al. [1^. Equation (38) of [l^ gives 
the cross section for semi-inclusive production of a pion 
0.5 1 1.5 0.5 1 1.5 2 
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• exp, Q'=2.3 GeV 
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FIG. 15: Ratio J^^^.^J^''^^^"'" for Kin2 (open circles) and Kin3 
(solid circles) plotted as a function of p,rx. Error bars rep- 
resent statistical errors only. We fitted to each kinematics a 
model by Anselmino et al. in [2^ using as a free 

parameter, where (fcx) is the intrinsic transverse momentum 
of quarks and {p\) is the transverse momentum transferred 
during the hadronization process. The reduced of the fits 
are 2.12 for Kin3 and 2.65 for Kin2. 



(valid for any hadron): 



they extracted these values in the inclusive region, im- 
plying a high multiplicity of particles, whereas in our 
data, the multiplicity of particles is unity. Typically, 
Anselmino et al. fit their model with data covering the 
range 0.1 < Zh < 1.0, with most of the statistics within 
Zfi < 0.4, whereas our data are within Zh > 0.9. Further- 
more, Kaskulov et al. [lB| used a value of 1.4 GeV^ for 
the rms transverse momentum of partons in their fit of 
the Hall C tt"*" data. The exclusive limit of SIDIS could be 
defined by a SIDIS-inspired model applicable to data at 
Ztt — >■ 1.0 or, more practically, when the measured hadron 
carries such a large fraction z^r of the total energy of the 
reaction that it does not allow the production of another 
particle. 

The HERMES and COMPASS collaborations have 
published (cos 2(/)^) moments of n'^ and 7r~ SIDIS, 
including up to 0.7 [sof . However, it is not possible to 
make a direct comparison to our axT tt" data as the 7r+ 
and 7r~ moments on the proton have different signs and 
magnitudes for Boer-Mulders effect. On the other hand, 
the higher twist Cahn effect, which also contributes to 
cttt, does not give by itself a satisfying description of 



dxBjdQ'^dz-^cPp^^ 



l + (l-y)2-4 



'-h{xB-Mz^) 



(2 - y)y/T^{k1)z^p^j 



cost 



1 



where y — pq/pk, and z^ = pp-n/pq is the fraction of 
the reaction energy carried by the measured hadron, 
and the quantities between angle brackets arc the stan- 
dard deviations of transverse momentum distributions, 
which are approximated as Gaussian, (fc^) stands for 
the parton transverse momentum in the proton, and 
^P^i) ~ (pi) + ^'ii^'i) is the measured transverse mo- 
mentum of the observed hadron, where (p^) stands for 
the transverse momentum of the hadron with respect to 
the direction of the struck quark. The idea is to adjust 
the ratio of cos over constant term in brackets of Eq. 



34l) by adjusting only the parameter (p\)/{k\ 



v/2£l(1 + e)(TTL 



4(2 - y)Vl - yzTrP7r± 



y? 



(35) 

Two conclusions arise from the fits shown in Figure 1151 
(1) the minus sign affecting the cos i^^ term in the SIDIS 
model is in agreement with the gtl and (2) {p\) must 
be equal to ^ 5.0 x (fcj ) to reproduce the data. 

The authors of [29| adjusted their model to semi- 
inclusive data. They give (/c^) = 0.25 GeV^ and {p\) = 
0.20 GcV^ giving a ratio {p\) / {k']) - 0.8. However, 



X. CONCLUSIONS 

We extracted the separated differential 7r° cross section 
at Jefferson Lab, Hall A, at four kinematic settings: Kin2 
and Kin3 with a 3% statistical precision, and KinX2 and 
KinXB with a 5% statistical precision. We studied the 
dependence of the hadronic tensor with the two first 
settings, and the xbj dependence with the latter two. 

The shape and order of magnitude of the cross section 
componants indicate that the t-channcl meson-exchange 
model is able to reproduce the total 7r° cross section, but 
it would still need improvement for the description of the 
other components. 

Table I VII summarizes the contradiction between our 
data and the leading twist QCD prediction for high . 
Instead of an ~ dependence we find, under the as- 
sumption that (Tt is negligible (which is very unlikely), 
& dependence for ctl. On the other hand, the cross 
section extracted may show an analogy with the formal- 
ism of SIDIS at the exclusive limit. Our ep — > epi:^ data, 
and the Hall C ep — >■ epir'^ data are important bases for 
studying the applicability of the SIDIS concepts to exclu- 
sive data. To improve the understanding of our data, we 
have run another tt'^ experiment in Fall 2010, at two beam 
energies, allowing us to disentangle Ci^fp from 
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dependence 


XBj dependence 




Kin3 Kin2 


KinX3 KinX2 


iVgcn 

fCdt 


15516 23429 
2.14 X 10^ 2.14 X 10^ 
5.10 X 10'-* nb"^ 2.99 x 10^ nb"^ 


5952 9860 


(GeV^) 

W (GeV) 
t^in (GeV") 
e 


2.350 ±0.002 1.941 ±0.010 
0.368 ± 0.001 0.368 ± 0.005 
2.217 ±0.004 2.055 ±0.012 
-0.173 ±0.001 -0.170 ±0.005 
0.649 ± 0.002 0.769 ± 0.003 


2.155 ±0.268 2.073 ±0.001 
0.335 ± 0.045 0.394 ± 0.003 
2.272 ±0.072 2.016 ± 0.008 
-0.137 ±0.048 -0.199 ±0.003 
0.648 ± 0.001 0.768 ± 0.003 


Eo (GeV) 
(GeV) 
q'"'' (GeV) 
p--- (GeV) 
(GeV) 


5.752 ±0.001 5.753 ±0.001 
2.348 ± 0.007 2.937 ± 0.020 
3.734 ± 0.007 3.143 ±0.017 
0.904 ± 0.002 0.806 ± 0.007 
0.910 ± 0.002 0.813 ±0.007 


5.752 ±0.001 5.753 ±0.001 
2.321 ±0.029 2.951 ±0.016 
3.732 ±0.009 3.151 ±0.014 
0.937 ±0.043 0.783 ±0.005 
0.942 ± 0.042 0.790 ± 0.005 



TABLE VII: Average quantities weighted with the cross section for the four kinematics of the experiment. Errors are the 
maximal deviation of the values in the seven fmin — t bins, compared to the averages hsted. 





dependence 




XBj dependence 


tmin-t (GeV") sin^^-'"' 


sin" e^ "" 


tmin-t (GeV") sin^^--"' 


sin" ei-""- 


Kin3 


KinA'3 


0.0095 


0.077 


0.007 


0.0095 


0.076 


0.007 


0.0298 


0.144 


0.021 


0.0297 


0.143 


0.020 


0.0546 


0.194 


0.038 


0.0545 


0.193 


0.037 


0.0844 


0.241 


0.058 


0.0843 


0.240 


0.058 


0.1188 


0.285 


0.081 


0.1188 


0.284 


0.081 


0.1583 


0.328 


0.108 


0.1579 


0.326 


0.106 


0.2063 


0.372 


0.139 


0.2057 


0.370 


0.137 


Kin2 


Kin A 2 


0.0094 


0.085 


0.008 


0.0094 


0.085 


0.008 


0.0296 


0.159 


0.026 


0.0296 


0.160 


0.026 


0.0541 


0.215 


0.046 


0.0542 


0.216 


0.047 


0.0839 


0.267 


0.071 


0.0840 


0.268 


0.072 


0.1179 


0.315 


0.099 


0.1181 


0.316 


0.100 


0.1576 


0.362 


0.131 


0.1579 


0.364 


0.133 


0.2050 


0.410 


0.168 


0.2051 


0.412 


0.170 



TABLE VIII: Values for tmin — t, sinO'^'^' and sin" e^ ™-, wei ghted by the cross section. 
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dependence 




Kin3 






Kin2 




SBj = 0.369, = 2.350GeV2 




0.368, Q2 = 1.941GeV^ 


fmin - t 


dar /dt + 


thdaL/dt 


GeV^ 




nb/GeV^ 




0.010 


377 ± 10 ± 12 






571 ± 10 ± 24 


0.030 


381 ± 12 ± 12 






600 ± 12 ± 25 


0.054 


403 ± 10 ± 13 






641 ± 12 ± 27 


0.084 


425 ± 11 ± 14 






673 ± 15 ± 28 


0.118 


418 ± 11 ± 14 






645 ± 16 ± 27 


0.158 


395 ± 13 ± 13 






636 ± 25 ± 27 


0.206 


384 ± 13 ± 13 






628 ± 36 ± 26 




duT h/dt 


0.010 


-13 ± 23 ± 10 






17 ± 19 ± 13 


0.030 


38 ± 26 ± 24 






-43 ± 22 ± 12 


0.054 


-25 ± 22 ± 11 






-23 ± 21 ± 12 


0.084 


-26 ± 25 ± 13 






-19 ± 27 ± 14 


0.118 


-75 ± 24 ± 9 






-103 ± 30 ± 21 


0.158 


-91 ± 30 ± 8 






-185 ± 52 ± 43 


0.206 


-123 ± 31 ± 10 






-189 ± 74 ± 34 




daTT 1 dt 


0.010 


-12 ± 23 ± 14 






-39 ± 19 ± 7 


0.030 


-25 ± 27 ± 15 






-110 ± 24 ± 13 


0.054 


-74 ± 22 ± 4 






-141 ± 22 ± 17 


0.084 


-64 ± 25 ± 14 






-174 ± 28 ± 17 


0.118 


-124 ± 24 ± 16 






-319 ± 29 ± 23 


0.158 


-137 ± 29 ± 15 






-352 ± 45 ± 53 


0.206 


-134 ± 30 ± 15 






-343 ± 57 ± 68 




doTL' /dt 


0.010 


9 ± 49 ± 20 






31 ± 51 ± 15 


0.030 


119 ± 55 ± 21 






136 ± 61 ± 24 


0.054 


129 ± 46 ± 12 






61 ± 56 ± 41 


0.084 


151 ± 51 ± 30 






123 ± 68 ± 20 


0.118 


153 ± 47 ± 17 






120 ± 69 ± 24 


0.158 


87 ± 54 ± 23 






142 ± 91 ± 36 


0.206 


127 ± 51 ± 15 






76 ± 99 ± 80 



TABLE IX: Separated cross-section values from Eq. (|19[l (first quoted value) with statistic errors (second quoted value) and 
systematic errors (third quoted value) for each of the seven considered bins. 
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XBj dependence 




KinX3 






KmX2 




XBj = 0.335, = 2.155GeV2 




0.394, Q2 = 2.073GeV^ 


fmin - t 


dar /dt + 


thdah/dt 


GeV^ 




nb/GeV^ 




0.010 


439 ± 19 ± 14 






635 ± 17 ± 26 


0.030 


437 ± 22 ± 14 






703 ± 21 ± 29 


0.054 


457 ± 18 ± 15 






683 ± 19 ± 28 


0.084 


442 ± 21 ± 14 






688 ± 23 ± 29 


0.118 


466 ± 22 ± 15 






682 ± 23 ± 28 


0.158 


407 ± 29 ± 13 






662 ± 34 ± 28 


0.205 


406 ± 34 ± 13 






591 ± 44 ± 25 




duT h/dt 


0.010 


20 ± 46 ± 38 






-26 ± 30 ± 22 


0.030 


2 ± 50 ± 17 






-100 ± 37 ± 61 


0.054 


-28 ± 43 ± 15 






-88 ± 32 ± 54 


0.084 


-37 ± 50 ± 19 






-68 ± 38 ± 487 


0.118 


-74 ± 55 ± 27 






-170 ± 40 ± 562 


0.158 


-188 ± 80 ± 27 






-155 ± 63 ± 657 


0.205 


-174 ± 90 ± 32 






-228 ± 82 ± 738 




daTT 1 dt 


0.010 


-16 ± 44 ± 16 






-63 ± 33 ± 18 


0.030 


-44 ± 50 ± 32 






-83 ± 41 ± 22 


0.054 


-63 ± 42 ± 15 






-153 ± 36 ± 24 


0.084 


-114 ± 47 ± 8 






-186 ± 43 ± 78 


0.118 


-156 ± 50 ± 18 






-327 ± 44 ± 109 


0.158 


-244 ± 66 ± 35 






-247 ± 65 ± 141 


0.205 


-124 ± 69 ± 42 






-444 ± 82 ± 183 




doTL' /dt 


0.010 


68 ± 97 ± 35 






-23 ± 84 ± 138 


0.030 


12 ± 109 ± 39 






112 ± 100 ± 104 


0.054 


236 ± 88 ± 19 






50 ± 90 ± 63 


0.084 


126 ± 99 ± 26 






211 ± 104 ± 95 


0.118 


119 ± 93 ± 22 






3 ± 106 ± 111 


0.158 


246 ± 106 ± 89 






78 ± 136 ± 126 


0.205 


177 ± 104 ± 30 






62 ± 146 ± 146 



TABLE X: Separated cross-section values from Eq. (|19p (first quoted value) with statistic errors (second quoted value) and 
systematic errors (third quoted value) for each of the first seven bins in fmin — t for 1.95 GeV^ < < 2.25 GeV^. 
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dependence 


Kin3 XBj = 


0.368, = 2.350 (GeV^) 




[-562 ± 62 ± 32] x sin^ d"^^' cos 2<1>^ nb 




[97 ± 18 ± 8] X sin B%^- cos $^ nb 


3m(W,,) = 


[-206 ± 35 ± 13] X sine^ "^- sin<l>^ nb 


Kin2 XBj = 


0.368, Q2 = 1.941 (GeV^) 


2 


[-1024 ± 58 ± 51] X sin^ e^ ™' cos 2$^ nb 




[82 ± 17 ± 11] X sin^^ '"- cos*^ nb 


9m(H/,,) = 


[-142 ±38 ±19] X sine^ ""- sin$^ nb 



TABLE XI: $.^-dependent hadronic tensor parametrization for constant xbj. The first error is the statistical error, the second 
is the systematic error. 



iBj dependence 



KinX3 XBj = 0.335, = 2.155 (GeV 



= [-770 ± 135 ± 63] X sin^ 6'^ ""- cos 2$^ nb 
3?e(W,,) = [121 ± 43 ± 17] X sinS^-"^' cos$^ nb 
'^miyV^^) = [-278 ± 69 ± 28] x sine*^'""' sin$^ nb 



KinX2 XBj = 0.394, = 2.073 (GeV^) 



— 2 - [-1003 ± 86 ± 153] X sin^ S^ ™' cos 2<1>^ nb 

Ke(W^^) = [163 ± 24 ± 72] x sin e*^ ""' cos $^ nb 
QmiW^^) = [-101 ± 58 ± 55] x sine*^-™' sin$^ nb 



TABLE XII: $.n--dependent hadronic tensor parametrization for constant Q^. The first error is the statistical error, the second 
is the systematic error. 
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